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ABSTRACT 


A new  method  has  been  developed  that  uses  satellite  along- 
track  position  errors  to  obtain  information  on  atmospheric  densi- 
ties at  altitudes  of  900  to  1200  km.  The  data  base  used  was  the 
dally  tracking  results  from  the  Navy  Navigation  Satellite  System 
(NNSS). 


Densities  derived  from  these  data  (for  the  years  1974-76) 

show  a pronounced  semiannual  variation,  with  the  April/October 
maxima  exceeding  the  August /January  minima  by  factors  of  2 to  4. 
The  variation  is  more  pronounced  than  was  previously  reported,  and 
double  the  best  current  modeled  variation.  The  current  resolution 
with  this  technique  is  estimated  to  be  1 x 10~19  gni/cm3,  or  about 
lOX  of  the  total  density  at  these  altitudes. 
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1.  SUMMARY 


We  have  devised  a method  that  employs  the  along-track  posi- 
tion errors  of  a satellite  system  to  correct  the  Jacchla  air- 
density  model.  We  have  used  the  Navy  Navigation  Satellite  System 
(NNSS,  or  Transit  System),  which  consists  of  a constellation  of 
satellites  in  near-circular  polar  orbits  at  altitudes  of  900  to 
1200  km.  The  satellite  along-track  error  is  very  sensitive  to 
drag  errors;  an  error  in  the  mean  density  of  2 x 10"T9  gm/cm^  will 
result  in  a day  in  an  along-track  error  that  will  grow  quadrati- 
cally  to  14  m.  This  is  easily  detectable  with  present  satellite- 
position  precision  (Ref.  1).  From  these  considerations  we  esti- 
mate that  the  resolution  with  the  technique  (using  1-day  aver- 
ages) is  about  1 X 10~19  gm/cm^  (10%  of  the  density  at  these 
altitudes) . 

Using  the  new  method,  we  analyzed  the  data  from  six  NNSS 
satellites.  Dally  values  of  satellite  along-track  position  errors 
were  converted  to  dally  mean-density  corrections,  which  were  then 
added  to  the  mean  modeled  density  (Ref.  2). 

The  semiannual  variation  in  the  upper  atmosphere  density  is 
clearly  evident  in  the  data.  Our  results  confirm  the  earlier 
Echo  2 (Refs.  3,  4,  and  5)  and  Calsphere-1  (Ref.  6)  findings  that 


Ref.  1.  H.  D.  Black,  R.  E.  Jenkins,  and  L.  L.  Pryor,  "The 
Transit  System,  1975,"  APL/JHU  TG  1305,  December  1976  (also  pre- 
sented at  the  56th  Annual  American  Geophysical  Union  Meeting, 

16-20  June  1975). 

Ref.  2.  L.  G.  Jacchia,  "Static  Diffusion  Models  of  the 
Upper  Atmosphere  with  Empirical  Temperature  Profiles,"  Smithsonian 
Contribution  to  Astrophysics.  Vol.  8,  No.  9,  1965. 

Ref.  3.  G.  E.  Cook  and  D.  W.  Scott,  "Exospheric  Densities 
Near  Solar  Minimum  Derived  from  the  Orbit  of  Echo-2,"  Planet . 

Space  Scl.  . Vol.  14,  pp.  1149-1165,  1966. 

Ref.  4.  G.  E.  Cook  and  D.  W.  Scott,  "Variations  in  Exo- 
spheric Density  at  Heights  Near  1100  km.  Derived  from  Satellite 
Orbits,"  Planet.  Space.  Scl..  Vol.  15,  pp.  1933-1956,  1967. 

Ref.  5.  G.  E.  Cook  and  D.  W.  Scott,  "The  Semi-Annual  Varia- 
tion on  Air  Density  at  Height  of  1100  km  from  1964  to  1967,"  Planet. 
Space  Sci. , Vol.  17,  pp.  107-119,  1969. 


Ref.  6.  G.  E.  Cook,  "The  Large  Semi-Annual  Variation  in  Exo- 
spheric Density:  A Possible  Explanation,"  Planet.  Space  Sci..  Vol. 
15,  pp.  627-632,  1967. 
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the  density  at  heights  of  approximately  1000  km  shows  a pronounced 
semiannual  variation.  The  1975  data  Indicate  an  October  maximum 
exceeding  the  July  minimum  by  a factor  of  2.5  and  an  April  maximum 
nearly  double  the  January  minimum.  Satellite  1973-81a  for  the 
years  1974-1976  shows  an  October  maximum  exceeding  the  July  mini- 
mum by  a factor  of  3.6.  The  discrepancy  may  be  accounted  for,  In 
part,  by  changing  the  modeled  diurnal  effect.  We  also  observed 
the  existence  of  local  July  maxima  In  the  data  from  satellites 
1973-81a  and  1970-67a.  The  maxima  may  be  partly  accounted  for  by 
an  error  In  the  amplitude  and/or  phase  of  the  diurnal  bulge  and/or 
the  seasonal-latitudinal  variations  of  helltim. 

Improvements  In  gravity  models  are  also  possible  with  satel- 
lites In  near  resonance  with  harmonics  In  the  geopotentlal  expan- 
sion. For  example,  the  data  from  satellite  1967-34a  clearly  ex- 
hibit a 28-day  resonant  oscillation  arising  from  errors  In  the 
27th-order  harmonics.  Future  NNSS  satellites  may  exhibit  other 
strong  resonances  as  well.  If  unexpected  resonance  effects  appear 
we  also  would  have  to  remove  these,  possibly  to  get  access  to  the 
density  variations. 
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2 . BACKGROUND 


During  the  period  1958-1963  we  worked  intensively  on  the 
analysis  and  software  for  the  Navy  Navigation  System  (Ref.  7). 

Our  effort  was  concerned  mostly  with  the  orbit  determination  pro- 
grams that  are  used  daily  in  computing  the  ephemerldes  of  the 
Transit  satellites.  The  five  satellites  are  all  in  polar  orbits 
at  different  nodal  longitudes  and  at  altitudes  of  950  to  1200  km. 

Since  1963,  when  the  program  became  operational,  APL  has 
been  responsible  for  updating  and  improving  the  precision  of  the 
system  via  software  changes.  Our  efforts  (of  geophysical  signifi- 
cance) have  concentrated  on  the  following:  (a)  Improvements  in 
the  geopotential  model  of  the  earth  Including  the  discovery  of 
"resonant"  geopotential  effects  (Refs.  8 and  9),  a revision  of  the 
Jacchia  air-density  model  to  adapt  it  to  the  above-800-km  altitude 
regime  (Ref.  10),  an  implementation  of  polar  motion  in  the  ephem- 
erls  construction  (Ref.  11),  and  construction  of  a model  of  the 
tropospheric  refraction  effect  (Refs.  7,  12,  and  13). 

The  data  we  have  used  in  our  air-density  studies  are  the 
normal  day-to-day  quality  control  statistic  on  the  ephemeris 


Ref.  7.  H.  D.  Black,  "Position  Determination  Using  the 
Transit  System,"  International  Geodetic  Symposium  on  Satellite 
Doppler  Positioning,  Las  Cruces,  NM,  12-14  October  1976. 

Ref.  8.  W.  H.  Guler,  "Geodetic  Problems  and  Satellite 
Orbits,"  Lectures  in  Applied  Mathematics,  Vol.  6,  Space  Mathe- 
matics, Part  II.  American  Mathematical  Society,  1966. 

Ref.  9.  S.  M.  Ylonoulls,  "Determination  of  Coefficients 
Associated  with  the  Geopotential  Harmonics  of  Order  Thirteen," 

J.  Geophys.  Res.,  Vol.  71,  No.  6,  March  1966,  p.  1768. 

Ref.  10.  A.  Eisner,  "Atmospheric  Density  Studies,"  APL/JHU 
TG  951,  December  1967. 

Ref.  11.  V.  L.  Plsacane,  B.  B.  Holland,  and  H.  D.  Black, 
"Recent  (1973)  Improvements  in  the  Navy  Navigation  Satellite  Sys- 
tem," Navigation.  Vol.  20,  No.  3,  Fall  1973,  pp.  224-229. 

Ref.  12.  H.  S.  Hopfleld,  "Two-Quartlc  Tropospheric  Refrac- 
tlvity  Profile  for  Correcting  Satellite  Data,"  J.  Geophys.  Res. , 
Vol.  74,  No.  18,  August  1969,  pp.  4487-4499. 

Ref.  13.  S.  M.  Ylonoulls,  "Algorithm  to  Compute  Tropospheric 
Refraction  Effects  on  Range  Measurements,"  J.  Geophys.  Res.,  Vol. 
75,  No.  36,  20  December  1970,  pp.  7636-7637. 
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production.  Prior  to  1967,  we  had  used  the  Harrls-Prlester 
model  of  the  upper  air  density  (Ref.  14)  In  computing  the  ephem- 
erldes.  The  rising  solar  cycle  forced  us  to  Implement  the  newer 
Jacchla  model  In  hopes  that  It  would  diminish  the  ephemerls 
error.  (The  normal  operation  of  Transit  requires  a prediction 
of  the  satellite  position  12  to  24  hours  Into  the  future.  Errors 
In  the  density  model  cause  period  errors  In  the  satellite  ephem- 
eris  that,  if  large  enough,  can  be  catastrophic.) 

We  Implemented  a strict  interpretation  of  the  Jacchla  alr- 
denslty  model  (Ref.  2)  In  the  ephemerls  complex,  realizing  that 
no  data  above  800  km  had  been  used  In  the  model  development. 

We  attempted  to  implement  it  in  the  dally  satellite  opera- 
tions but  it  failed.  We  later  learned  that  we  had  attempted  the 
implementation  during  a very  severe  magnetic  storm  (May  1967).  An 
intensive  investigation  (Ref.  10)  showed  that  the  model  consis- 
tently overestimated  the  air  density  during  disturbed  magnetic 
periods.  Minor  parameter  adjustments  to  the  basically  sound 
Jacchla  model  solved  the  problem.  We  have  used  It  ever  since 
with  dally  Insertion  of  solar  and  magnetic  Indexes. 

About  18  months  ago,  we  shifted  the  Transit  system  from  the 
APL  4.5  geopotential  model  (Ref.  15)  to  the  WGS-72  model  (Ref.  16). 
For  some  time  prior  to  this  we  had  noticed  an  apparent  anomaly  in 
the  results  from  one  of  the  satellites.  Of  the  five  satellites, 
the  anomalous-behaving  satellite  had  the  lowest  perigee  height. 

We  had  sought  In  vain  for  the  cause  of  the  anomaly.  When  the 
anomaly  persisted  after  the  change  In  geopotential  models,  we  seri- 
ously suspected  the  drag  model.  One  of  us  (Eisner)  remembered  that 
the  semiannual  term  in  the  Jacchla  model  phase-matched  our  data; 
this  was  the  key  to  resolving  the  problem. 

The  Jacchla  model  (Ref.  2)  Incorporated  a semiannual  varia- 
tion in  the  mean  exospheric  temperature  with  an  amplitude  propor- 
tional to  the  average  10.7-cm  solar  radiation  flux,  Fio.7»  Cook 
(Ref.  17)  found  that  the  Jacchla  model  under-represents  the 


Ref.  14.  I.  Harris  and  W.  Prlester,  "Theoretical  Models  for 
the  Solar-Cycle  Variation  of  the  Upper  Atmosphere,"  J.  Geophys.  Res., 
Vol.  67,  No.  12,  November  1962,  pp.  4585-4591. 

Ref.  15.  H.  D.  Black,  "Doppler  Tracking  of  Near-Earth  Satel- 
lites," APL/JHU  TG  1031,  1968. 

Ref.  16.  T.  0.  Seppelln,  "The  Department  of  Defense  World 
Geodetic  System  1972,"  The  Canadian  Surveyor.  Vol.  28,  No.  5, 

Ottawa,  Canada,  December  1974,  pp.  496-506. 

Ref.  17.  G.  E.  Cook,  "The  Semi-Annual  Variation  In  the  Upper 
Atmosphere:  A Review,"  Ann,  de  Geophys.,  Vol.  25,  1969,  pp.  451-469. 
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semiannual  effect  at  heights  above  1000  km,  a conclusion  that  was 
confirmed  in  our  investigations  of  the  NNSS  satellites.  Jacchia 
(Ref.  18)  acknowledged  the  shortcomings  of  the  earlier  model  and 
switched  from  a solar-actlvlty-dependent  effect  (caused  by  tempera- 
ture variations)  to  a functional  form  that  is  strongly  height- 
dependent  but  is  Independent  of  solar  activity.  In  his  latest 
model  he  expresses  the  semiannual  density  variation  in  the  form 


A logj^Q (semiannual)  “ f(z)  g(t) 


(1) 


where  f(z)  relates  the  amplitude  to  height  and  g(t)  represents  the 
average  density  variations  as  a function  of  time. 

Helium  is  the  dominant  component  at  1000  km.  Hedln  et  al. 
(Ref.  19)  found  a morning  maximum  for  helium  located  in  the  winter 
hemisphere.  Von  Zahn  et  al.  (Ref.  20)  emphasize  the  dominance  of 
the  annual  migration  of  the  helium  bulge  over  the  diurnal  variation 
and  the  difficulty  of  separating  the  two  effects.  The  so-called 
"winter  helium  bulge,"  first  discovered  in  1967  (Ref.  21),  was  later 
shown  to  be  a special  case  of  the  more  general  seasonal-latitudinal 
variation  in  density.  Mayr  and  Volland  (Ref.  22)  concluded  from 
all  the  available  evidence  that  the  semiannual  effect  consists  of 
two  components  (a  global  component  and  a latitude-dependent  compo- 
nent), both  of  which  are  associated  with  temperature  variations. 


Ref.  18.  L.  G.  Jacchia,  "Revised  Static  Models  of  the 
Thermosphere  and  Exosphere  with  Empirical  Temperature  Profiles," 
Smithsonian  Astrophyslcal  Observatory  Special  Report  332,  May  1971. 

Ref.  19.  A.  E.  Hedin,  H.  G.  Mayr,  C.  A.  Reber,  and  N.  W. 
Spencer,  "Empirical  Model  of  Global  Thermospheric  Temperature  and 
Composition  Based  on  Data  from  OGO-6 , Quadrapole  Mass  Spectrometer," 
J.  Geophys.  Res.,  Vol.  70,  No.  1,  1975,  pp.  215-225. 

Ref.  20.  V.  Von  Zahn,  W.  Kbehnleon,  K.  H.  Frlcke,  V.  Laux, 

H.  Trlnds,  and  H.  Volland,  "ESRO  4 Model  of  Global  Thermospheric 
Composition  and  Temperatures  During  Times  of  Low  Solar  Activity," 
Geophys.  Res.  Lett..  Vol.  4,  No.  1,  January  1977,  pp.  34-36. 

Ref.  21.  G.  M.  Keating  and  E.  J.  Prior,  "The  Winter  Helium 
Bulge,"  Space  Res.,  Vol.  8,  1968,  pp.  982-992. 

Ref.  22.  H.  G.  Mayr  and  H.  Volland,  "Theoretical  Model  for 
the  Latitude  Dependence  of  the  Thermospheric  Annual  and  Semi-Annual 
Variations,"  J.  Geophys.  Res.,  Vol.  77,  No.  34,  1972,  pp.  6774-6790. 
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Jacchla's  1965  model  was  heavily  weighted  toward  the  lower 
regions  of  the  atmosphere  since  few,  If  any,  data  were  available 
above  1000  km.  Jacchla's  latest  model  (Ref.  18)  has  Incorporated 
data  from  the  Echo-2  and  Calsphere-l  satellites  together  with 
seasonal-latitudinal  helium  variations  (Ref.  21),  thereby  extending 
the  model  well  Into  the  altitude  regime  of  the  Transit  satellites. 


12 


TH£  JOHNS  HOPKINS  UNlVt^SlTV 

applied  physics  LABOMATORY 


I 


4- 


1 

! 

I 


3.  CURRENT  STATUS 


SATELLITE  DATA  DESCRIPTION 

Each  satellite  In  the  constellation  of  polar  satellites 
transmits  Its  current  position  In  addition  to  the  doppler  carriers 
(a  frequency  pair  at  150  and  400  MHz) . The  position  is  a real-time 
readout  from  a predicted  ephemerls  that  Is  contained  In  the  satel- 
lite memory.  The  predicted  ephemerls  is  generated  as  follows.  On 
a continuing  basis,  a network  of  ground  stations  takes  all  possible 
doppler  data  from  passes  (transits)  of  the  satellites.  The  data 
are  then  transmitted  to  a central  processing  center.  There  they 
are  accumulated  into  36-hour  spans  for  each  satellite.  A least- 
squares  estimation  is  performed  on  each  span  to  determine  a set  of 
satellite  initial  conditions  that,  when  used  to  generate  the  satel- 
lite ephemerls,  gives  the  best  fit  to  the  doppler  data.  (Numerical 
integration  of  the  forces  acting  on  the  satellite  is  used  for  this 
computation.)  The  set  of  initial  conditions  is  then  used  to  pre- 
dict the  satellite  ephemerls  approximately  30  hours  beyond  the  end 
of  the  data  span.  The  predicted  ephemerls  is  then  transmitted  to 
the  satellite,  where  it  is  stored  in  the  satellite  memory. 

The  input  or  "starting"  orbit  for  a current  36-hour  span  is 
obtained  by  re-epochlng  the  orbit  determined  from  the  previous  36- 
hour  span.  An  analysis  of  the  doppler  residuals  over  the  new  span 
measures  the  quality  of  the  predicted  ephemerls  currently  being 
used  in  the  satellite.  The  rms  of  the  residuals,  T (in  meters),  is 
plotted  versus  time  in  Figs.  1,  2,  and  3.  The  changes  in  initial 
conditions  obtained  in  the  estimation  procedure  also  reflect  how 
the  accuracy  of  the  predicted  ephemerls  degrades  with  time.  One  of 
the  initial-condition  parameters  is  highly  sensitive  to  secular  or 
long-period  errors  in  the  satellite-motion  model.  The  parameter 
is  routinely  plotted  to  provide  a "quick-look"  estimate  of  the  pre- 
dlcted-ephemerls  precision.  It  provides  an  excellent  data  source 
for  studying  the  long-term  behavior  of  atmospheric  density  model 
errors.  This  is  shown  in  Appendix  A. 

From  Eq.  A-12  of  Appendix  A we  obtain 
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Day  number  j _ rms  satellite  position  error 

Ai  — rate  of  growth  of  satellite 
along-track  error 


Fig.  1 NNSS  Satellite  Summary  Statistics,  Satellite  1973-81a  (30200) 
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where 

6p  ■ error  in  atmospheric  density 

n ■ satellite  mean  motion 

■ coefficient  of  drag 
A 

— - area-to-mass  ratio  of  the  satellite 

V ■ satellite  speed  relative  to  the  atmosphere 

tg.t  ■ starting  epochs  associated  with  consecutive 

fitting  Intervals 

Equation  A-12  Is  the  formula  used  to  convert  the  (rate  of 
growth  of  satellite  along-track  error)  data  Into  mean  density 
measurements . 

The  A]^  parameter  can  also  be  expressed  as  a change  In  the 
satellite's  semi-major  axis,  a,  at  t ■ t]^  via 


(see  Figs.  1,  2,  and  3). 


DISCUSSION  OF  INITIAL  RESULTS 

Our  attention  was  Initially  drawn  to  the  dally  summary 
statistics  for  NNSS  satellite  1973-81a  (Fig.  1).  Peak  errors  oc- 
curring on  Day  95  (6  April)  and  Day  295  (23  October)  brought  to 
mind  the  semiannual  variation  In  the  atmosphere.  Summary  statis- 
tics for  the  other  five  satellites  In  the  NNSS  (Table  1)  reinforced 
our  Initial  Impression  that  we  were  dealing  with  atmospheric  semi- 
annual variation  model  errors  In  the  data  (Figs.  2 and  3).  Satel- 
lite 1967-34a  (Fig.  3)  exhibits  the  effects  of  errors  In  the 
modeled  semiannual  variation  that  are  somewhat  corrupted  by  an  os- 
cillation with  a period  of  about  28  days.  (This  particular  satel- 
lite Is  strongly  resonant  with  the  27th-order  harmonics  In  the 
geopotentlal;  In  Fig.  3 we  see  the  residual  errors  In  these  terms. 

We  have  developed  an  algorithm  for  converting  the  raw  data, 
A]^,  Into  an  equivalent  density  correction  to  the  modeled  mean  den- 
sity. We  replaced  the  dally  values  of  A^^  (which  are  somewhat  noisy, 
as  shown  In  Figs.  1,  2,  and  3)  with  10-day  averaged  values  that 
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Years 

Covered 

1976 

1976 

1976 

1976 

1976 

1974-1976 

Perigee 

Precession 

Rate 

(deg/day) 

-2.90 

-2.86 

-2.87 

-2.87 

-2.85 

-2.97 

Nodal 

Precession 

Rate 

(deg/day) 

0.021 

-0.038 

-0.077 

-0.002 

-0.009 

-0.019 

Perigee 

Position 

1 Jan  76 
(deg) 

180 

133 

119 

294 

37 

131 

Nodal 

Position 

1 Jan  76 
(deg) 

23 

- 47 

- 35 

- 84 
-109 

123 

Inclination 

(deg) 

90.20 

89.60 

89.25 

89.98 

90.10 

90.16 

Eccentricity 

1 

0.002 

0.002 

0.005 

0.008 

0.018 

0.017 

Semi- 

Major 

Axis 

(km) 

7449 

7463 

7454 

7462 

7464 

7399 

Satellite 

1 

1967-34a 

7 1967  48a 

1967  92a 

1968 12a 

1970-67a 

1973  81a 
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are  well-suited  to  the  study  of  a variation  having  a 180-day  period. 
The  resulting  density  corrections  for  satellite  1973-81a  for  the 
years  1974-1976  and  for  the  other  five  NNSS  satellites  for  the 
year  1976  are  presented  in  Figs.  4 and  5.  We  then  added  the  cor- 
rections to  the  modeled  mean  densities  (used  in  the  generation  of 
the  ephemerls).  Figures  6 and  7 present  the  actual  mean  densities 
experienced  by  the  six  NNSS  satellites.  Satellite  1973-81a  ex- 
hibits a pronounced  semiannual  variation  throughout  the  period 
studied  (1974-1976).  There  is  generally  good  agreement  in  both 
phase  and  amplitude  from  year  to  year. 

We  note  a curious  peaking  in  density  in  July,  particularly 
in  the  1976  data.  The  peaking  is  partially  accounted  for  by  the 
diurnal  bulge  "beating"  against  perigee  position:  The  NNSS  satel- 
lites have  nearly  polar  orbits;  consequently  their  orbital  planes 
precess  very  slowly  in  inertial  space  (Table  1) . Since  the  sun 
and  bulge  move  about  l®/day,  it  follows  that  the  satellite  plane 
contains  the  diurnal  bulge  twice  yearly  (once  each  on  the  ascend- 
ing and  descending  nodes  of  the  orbit).  The  net  effect  is  very 
smal.l  unless  the  eccentricity  is  large  enough  for  the  position  of 
perigee  relative  to  the  diurnal  bulge  to  play  a part.  The  diurnal 
and  semiannual  variations  for  1973-81a  and  1970-67a  (the  two  rela- 
tively eccentric  NNSS  satellites)  interfere  constructively  (or 
destructively)  resulting  in  the  appearance  of  an  enhanced  (or  sup- 
pressed) pseudosemiannual  variation.  Satellite  1973-81a  is  particu- 
larly interesting  (Fig.  8).  At  the  start  of  1976,  perigee  (located 
on  the  ascending-node  side  of  the  orbit)  moves  through  the  diurnal 
bulge,  resulting  in  peak  mean  densities.  One  hundred  eighty  days 
later,  perigee  has  precessed  540“  and  is  once  again  traversing 
through  the  diurnal  bulge  (on  the  descending-node  side  of  the 
orbit).  The  cycle  is  completed  180  days  later  (at  the  end  of  1976) 
when  the  situation  is  the  same  as  it  was  at  the  start  of  1976.  This 
pattern  is  repeated  over  several  years  because  of  the  very  small 
nodal  precession  rate  of  1973-81a  and  the  fact  that  perigee  returns 
to  its  starting  position  every  365  days.  The  peaking  densities  at 
the  beginning,  middle,  and  end  of  the  year  due  to  the  diurnal  bulge 
destructively  interfere  with  the  semiannual  minima  occurring  approxi- 
mately at  the  same  times.  The  result  is  particularly  evident  in 
July  and  August  where  an  expected  minimum  (particularly  in  1976 
data)  is  replaced  with  a local  maximum.  We  will  discuss  this  in 
more  detail  later. 

In  addition  to  the  semiannual  and  diurnal  variations,  we 
also  have  to  contend  with  the  seasonal  variations  in  the  concentra- 
tion of  helium  (Refs.  20  and  21).  This  effect  results  in  peak  den- 
sities occurring  in  January  and  July  (Fig.  8) . These  maxima  de- 
structively interfere  with  the  semiannual  minima  occurring  at  about 
the  same  time,  resulting  in  the  appearance  of  a reduced  semiannual 
amplitude. 
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p = 898  km  Day  number 

a = 1145  km 

Fig.  4 Corrections  to  Modeled  Atmospheric  Density  versus  Time  (derived  from 
Satellite  1973-81a  (30200)  data) 
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Fig.  8 Mean  Modeled  Densities  versus  Time  (modeled  semiannual  variation 
removed;  various  satellites) 
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Figures  6 and  7 clearly  illustrate  a needed  change  to  the 
1965  semiannual  variation  Jacchla  model  (Ref.  2).  This  need  had 
been  noted  (Refs.  3 and  17)  and  has  resulted  in  a functional  form 
Independent  of  solar  activity  (Ref.  18).  We  employed  this  model 
to  generate  densities  for  the  six  satellites.  Figures  9 and  10 
present  the  results.  Although  the  agreement  with  the  data  is 
much  better,  it  is  immediately  apparent  that  the  amplitude  of  the 
modeled  semiannual  variation  is  not  sufficiently  large. 

Our  next  step  was  to  estimate  the  semiannual  variation  from 
our  data.  We  removed  all  other  variations  from  our  data  using  the 
1965  model  (with  1971  diurnal  parameters)  but  omitting  the  modeled 
semiannual  variation  (Fig.  8).  The  results  are  plotted  in  Figs.  11 
and  12  and  summarized  in  Table  2.  We  note  a number  of  interesting 
features  in  these  plots: 

1.  A mid-year  flattening  and  possible  peaking  in  both  1973- 
81a  (all  3 years)  and  1970-67a  data, 

2.  A 28-day  oscillation  superimposed  on  the  semiannual 
variation  in  the  1967-34a  data,  and 

3.  A suggestion  of  peak  densities  in  May  (satellites  1967- 
48a  and  1968-12a) . 

The  first  of  the  above  (pronounced  in  the  1976,  1973-81a 
(30200)  data)  is  Interesting  and  puzzling.  Simply  raising  the  am- 
plitude of  the  diurnal  effect  Improves  the  situation  in  mid-year 
but  at  the  expense  of  a poorer  fit  of  the  data  at  the  beginning 
and  end  of  the  year.  The  unaccounted-for  seasonal  migration  of 
the  helium  bulge  likewise  will  tend  to  counteract  the  falling  den- 
sities in  January  and  July  and  give  the  appearance  of  flattened 
semiannual  minima.  An  interesting  possibility  that  comes  to  mind 
is  the  case  where  the  semiannual  effect  was  a result  of  the  varia- 
tion in  exospheric  temperature  (Ref.  2).  Under  low  solar  activity, 
minimum  nighttime  temperatures  during  the  July  semiannual  minimum 
may  fall  below  700  K,  a region  where,  for  altitudes  of  1000  km, 
density  actually  increases  with  falling  temperatures  (Fig.  13). 

The  same  local  July  minima  were  also  evident  in  the  Echo-2  derived 
densities  for  1965  and  1966  (Fig.  11  of  Ref.  4). 

The  second  point  was  mentioned  earlier  and  is  a resonance 
effect  (Ref.  23)  as  a result  of  errors  in  the  27th-order  modeled 
geopotential  coefficient.  This  once  again  points  out  the  high 


Ref.  2.  S.  M.  Ylonoulls,  "A  Study  of  the  Resonance  Effects 
Due  to  the  Earth's  Potential  Function,"  J.  Geophys.  Res.,  Vol.  70, 
No.  24,  December  1965,  pp.  5991-5996,  and  Vol.  71,  No.  4,  January 
1966,  pp.  1289-1291. 


Fig.  9 Mean  Atmospheric  Density  versus  Time  (modified  model;  from 
Satellite  1973-81a  (30200)  data) 


TMf  JOHNS  M0»*HINS  university 

APPLIED  PHYSICS  LABORATORY 

lAOHCl  MAHviAND 


THC  JOHNS  HOPKINS  UNIVERSITY 

APPLIED  PHYSICS  LABORATORY 

IAU#FI  MAn>'(AM0 

P 

Po 

— Mean  density  derived  from  NNSS 

- Modeled  mean  density  with  semiannual 
variation  removed 

2 

-1*  1 

^ ^ ^ 

la. 

hp  = 953  km 
ha  = 1221  km 

Satellite  1970  67a  (30190) 

, = 1030  km 
= 1140  km 


Satellite  1968- 12a  (30180) 


1040  km 
1115km 


Satellite  1967-92a  (30140) 


1071  km 
1101  km 


Satellite  1967-48a  (30130) 


1048  km 
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Satellite  1967-34a  (30120) 

I I I 


Day  (1976) 

Fig.  12  Semiannual  Variation  in  Density  versus  Time  (derived  from  data 
from  various  satellites) 
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Table  2 

Phase  and  Amplitude  of  the  Semiannual 
Variation  Derived  from  NNSS  Data 


Satellite 

Year 

Position  of  Maximum/Minimun 
(day  number) 

Density  maximum 
Density  minimum 

January 

(min) 

April 

(max) 

July 

(min) 

October 

(max) 

April 

(max) 

October 

(max) 

1967  34a 

25 

90 

190 

310 

2.03 

2.31 

1967  48a 

60 

95 

200 

285 

1.41 

2.20 

1967  92a 

1976 

30 

100 

195 

295 

1.60 

4.20 

1968  12a 

15 

65 

220 

280 

2.10 

1.75 

197a67a 

20 

85 

190 

290 

1.40 

1.60 

1974 

1 

80 

205 

290 

2.60 

4.80 

1973-81a 

1975 

1 

90 

205 

290 

2.80 

5.70 

1976 

1 

90 

180 

295 

2.46 

3.16 

Six-satellite 

Jan 

Mar 

July 

Oct 

mean  for  1976 

24 

28 

14 

19 

1.83 

2.54 

Satellite  1973-81a 

Jan 

Mar 

July 

Oct 

2.62  1 

4.55 

mean  for  1974-76 

1 

1 

27 

1 

14 

17 

sensitivity  of  this  type  of  data  for  detecting  effects  with  along- 
track  amplitude  of  10  to  IS  m and  periods  as  short  as  1 month. 

We  have  no  explanation  for  point  3 above.  Both  satellites 
are  nearly  circular  so  that  the  diurnal  effect  Is  of  little  conse- 
quence. We  cannot  Ignore  the  dip  In  Day  120  in  the  1967-48a  data 
and  on  Day  110  In  the  1968-12a  data  as  simply  noise  because  It  Is 
clearly  present  In  the  raw  dally  data. 
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Appendix  A 


DERIVATION  OF  Aj^-6p  RELATIONSHIP 


Using  a first-order  perturbation  theory  and  Ignoring  eccen- 
tricity effects  (Ref.  8),  the  equations  of  motion  governing  the 
propagation  of  satellite  orbit  errors  can  be  written  as 


H - 

3nH  - 2nL  - F„ 
rl 

altitude 

(A-la) 

L - 

2nH  - Fj^ 

along-track 

(A-lb) 

C + 

cross-track, 

(A-lc) 

where  H,  L,  and  C are  components  of  the  satellite  position  error 
resolved  (a)  In  the  direction  of  the  satellite  radius  vector, 

(b)  In  the  direction  of  the  satellite  velocity  (along-track) , and 

(c)  In  the  direction  of  the  angular  momentrum  vector  (cross-track) , 
respectively.  Errors  in  the  forces  acting  on  the  satellite  In 
these  three  directions  are  given  by  F^,  Fl,  and  F(^.  The  symbol  n 
Is  the  mean  motion  of  the  satellite.  Dots  denote  derivatives  with 
respect  to  time. 

The  current  knowledge  of  the  geopotential  dictates  that  the 
major  error  source  In  predicting  satellite  orbits  Is  In  the  drag 
force  model.  A mean  error  In  the  modeled  drag  force  Is  equivalent 
to  a constant  force  acting  In  the  along-track  direction  of  the 
satellite's  motion  (l.e.,  F^  “ SF^j  « constant).  The  motion  to 
Eq.  A-1  for  this  constant  force  is  given  by 


j cos  n(t  - tjj) 


H ■ j j sin  n(t  - t^) 


(t  - tg) 


L - Bq  + B^n(t  - tg)  + B^  cos  n(t  - t^)  + B^  sin  n(t 

+ B3  sin  n(t  - tg)  - I F^(t  - tp)^ 

C - B^  cos  n(t  - tg)  + Bj  sin  n(t  - t^)  . 


(A-2a) 


(A-2b) 


(A-2c) 
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The  Bj^(l“0, . . .5)  are  the  constants  of  Integration  and  reflect  the 
effects  of  Initial  condition  errors  at  t ■ tQ* 

From  Eq.  A-2  we  see  that  the  major  effect  of  a mean  error 
in  the  atmospheric  drag-force  model  is  to  cause  a quadratic  growth 
in  the  satellite's  along-track  error.  It  is  also  clear  that  in 
determining  corrections  to  the  initial  condition  parameters,  Bq 
and  Bj^  will  be  the  most  sensitive  to  the  Fj  error.  The  remaining 
parameters  are  more  sensitive  to  orbital  frequency  errors. 

To  estimate  the  contribution  of  6F(i  to  the  Bq  and  Bj^  parame- 
ters we  can  simplify  the  problem  by  (a)  minimizing  only  the  L 
residuals  and  (b)  using  only  the  Bq  and  Bi  parameters  in  the  theo- 
retical model.  Thus  we  minimize  the  function 


F - 


T E 2 
(L  - L ) dt 


(A- 3) 


with  respect  to  Bg  and  Bj^  where 


J - Bq  + B^  n(t  - tg)  (A-4) 

- - I «F^(t  - tg)^  (A-5) 


(cf  Eq.  A-2b)  and  tg  and  t2  are  the  beginning  and  ending  times, 
respectively,  associated  with  the  fitting  (or  "tracking")  span. 

The  least-squares  solution  for  Bg  and  Bj^  is  found  to  be 


’o  ■ i <'2  • ’’d 


“l  ■ ■ 2 *'2  ' '0*  n 


(A-6a) 


(A-6b) 
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Thus  the  residuals,  after  tracking,  are  given  by 


L 


R 


(t. 


tl)(t  - tj^,)  + (t  - 


(A-7) 


and  are  shown  In  Fig.  A-1.  The  next  fitting  Interval  covers  the 
time  span  from  ti  to  t3  and  the  residuals  to  be  minimized  are  given 
by  Eq.  A-7.  The  theoretical  model  used  to  fit  to  these  residuals 
Is  again  given  by 


+ Aj  n(t  - 


^1^ 


(A-8) 


where  the  epoch  is  now  at  tj^.  The  solution  obtained  for  Aq  and  A^ 
Is 


(A-9a) 


(A-9b) 


Equation  A-9b  gives  the  desired  dependence  of  the  Aj  Initial-condi- 
tion parameter  with  respect  to  a mean  drag-force  error. 

The  atmospheric  drag-force  acting  on  the  satellite  is 
modeled  as 


1 A „2 


(A-10) 


where 


A 

m 

V 

P 


coefficient  of  drag, 

area-to-mass  ratio  of  the  satellite, 

satellite  velocity  relative  to  the  atmosphere,  and 

atmospheric  density. 
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Satellite  along  track  error. 
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The  atmospheric  density  Is  computed  using  a modified  Jacchla  model 
(Ref.  10).  From  Eq.  A- 10  we  then  relate  6Fj  to  a mean  error  In 
the  computed  density,  5p: 


F 


d 


1 

2 


C - 
d m 


6p  . 


(A-11) 


Substituting  Eq.  A-11  into  Eq.  A-9b  and  solving  for  6p  yields 


6p 


n 


S ''"'4  - 'o> 


(A-12) 


Equation  A-12  Is  the  formula  used  to  convert  the  A^^  data  Into 
mean  density  measurements. 

The  A]^  parameter  can  also  be  expressed  as  a change  In  the 
satellite's  semi -major  axis,  a,  at  t > tj.  via 


A^  - - I 6a  . (A-13) 


Aq  - Fitted  along-track  error  at  t\ 
Ai  — Slope  of  line  fitted  to  L 


Fig.  A-1  Satellite  Along-Track  Errors  Induced  by  Drag 
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GLOSSARY 


Along-Track  Error 

The  difference  between  an  a priori  computed  and  a least- 
squares  estimated  satellite  position  Is  resolved  In  three  orthog- 
onal directions.  One  of  these  directions  lies  in  the  direction  of 
the  satellite  velocity  vector  at  a particular  Instant.  The  corre- 
sponding position  error  component  Is  called  "along-track. " 

APL  4.5  Geopotential  Model 

A representation  of  the  geopotential  (gravity  field)  as  a 
16th  order/degree  Legendre  polynomial  expansion.  This  model  was 
developed  at  APL  In  1966  and  has  been  made  obsolete  by  more  accu- 
rate representation. 

NNSS  — Navy  Navigation  Satellite  System 

Frequently  called  the  Transit  System.  The  system  consists 
of  five  polar  orbiting  satellites,  all  roughly  at  1000-km  altitude, 
a network  of  four  tracking  sites  (all  In  the  USA),  a computation 
center  In  California,  and  an  injection  station  to  transmit  ephem- 
eris  data  to  the  satellite. 

Resonant  Geopotential  Effects 

If  the  satellite  nodal  period  Is  (close  to)  an  integral  sub- 
multiple of  the  sidereal  day  then  small  errors  In  a corresponding 
sectorial  harmonic  can  cause  large  satellite  position  errors. 

WGS-72 


The  World  Geodetic  System  1972  is  a representation  of  the 
geopotential  developed  by  the  Department  of  Defense  In  1972 
(Ref.  16). 
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